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ABSTRACT 

We present very deep Wide Field Camera 3 (WFC3) photometry of a massive, compact galaxy 
located in the Hubble Ultra Deep Field. This quiescent galaxy has a spectroscopi c redshift z = 
1.91 and has been identified as an extremely compact galaxy by iDaddi et al.l (|2005). We use new 
f/pi60W imaging data obtained with Hubble Space Telescope/ WFC3 to measure the deconvolved 
surface brightness profile to H ss 28 mag arcsec -2 . We find that the surface brightness profile is well 
approximated by an n = 3.7 Sersic profile. Our deconvolved profile is constructed by a new technique 
which corrects the best-fit Sersic profile with the residual of the fit to the observed image. This allows 
for galaxy profiles which deviate from a Sersic profile. We determine the effective radius of this galaxy: 
r e = 0.42 ± 0.14 kpc in the observed i?pi6ow-band. We show that this result is robust to deviations 
from the Sersic model used in the fit. We test the sensitivity of our analysis to faint "wings" in the 
profile using simulated galaxy images consisting of a bright compact component and a faint extended 
component. We find that due to the combination of the WFC3 imaging depth and our method's 
sensitivity to extended faint emission we can accurately trace the intrinsic surface brightness profile, 
and that we can therefore confidently rule out the existence of a faint extended envelope around the 
observed galaxy down to our surface brightness limit. These results confirm that the galaxy lies a 
factor ~ 10 off from the local mass-size relation. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

A significant fraction of massive galaxies at 2 w 
2 are early-type galaxies containing quiescent stellar 
populations (e.g., Franx et al.l 120031 : IDaddi et"all 120051 : 
iKriek et al.l 120061 ) . These galaxies must have formed 
very early in the universe's history and can therefore 
provide important constraints on galaxy formation and 
evolution models. Many of these quiescent galaxies 
have been found to be extremely compact, with effec- 
tive radii a fac tor ~ 6 smaller than their low-z coun- 
terparts (e.g.. IDaddi et al.l [20051 : iTruiillo et all [2001 
van Dokkum et al.l 120081) . This is quite puzzling, since 
these compact galaxies are passively evolving and are 
therefore not expected to change st rongly in size or mass 
if they do not merge. We note that lMancini et al.l (|2010l ) 
find some large massive quiescent galaxies at z ~ 1.5, 
showing that not all massive quiescent galaxies at high 
redshift are compact. 

Within the context of current-day models, galaxy 
merg ers play an importan t role in galaxy evolution 
fe.g.. IWhite fc Frenklll99l . These mergers may cause 
compact z ~ 2 galaxies to grow "inside-out", i.e., 
the mergers would increase the size of the galaxies 
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(e.g., Ivan Dokkum et al.l 120101: IHopkins et all I2009bl) . 

Whether the resulting size growth is large en ough, how- 
ever, is uncertain (e.g.. iBezanson et al.ll2009D . 

Several authors have emphasised that there are sev- 
eral systematic uncertainties that affect both radius and 
mass determinations. First, effective radii may be un- 
derestimated due to complex morphologies. Specifically, 
an extended low surface brightness component could re- 
main undetected due to low signal-to- noi se ratio (S/N) 
thereby lowering the observ ed size ( e.g.. IHopkins et ah 
2009at iMancini et al.l 120101 bu t see [van Dokk um et al.l 
20081: Ivan der Wei et al.l 12008ft . Second, mass-to-light 
gradients may result in a luminosity-weighted effective 
radius that is smaller than the mass-weighted effectiv e 
radius (e.g.. IHopkins et al.l[2009al : IHopkins et al.ll2009bf 1. 
Such gradients arise in certain models for the forma- 
tion of massive e llipticals (e.g., iRobertson et al.l 120061 : 
iNaab et al.l 120071 ) . Finally, the inferred stellar masses 
may be affected by incorrect assumptions regarding the 
initial mass function (IM F) and stellar evolution models 
(e.g.. IMuzzin et aLll2009l and references therein). 

In this Letter we use new very deep near-infrared (NIR) 
imaging data from the Hubble Space Telescope's Wide 
Field Camera 3 (HST/WFC3) to investigate the possi- 
bility of size underestimation due to the lack of S/N. 
We examine the possibility of a "hidden" faint extended 
component being present in z ~ 2 compact quiescent 
galaxies, focusing on the most massive quiescent galaxy 
in the Hubble Ultra Deep Field (HUDF) (B eckwith et al 
2006) , which has previously been studied by IDaddi et al.l 
(2005). We adopt the following values for cosmological 
parameters: H = 70 km s _1 Mpc -1 , Q m — 0.3, and 
Q\ = 0.7. All stellar masses are derived assuming a 
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Kroupa IMF (jKroupal [200lh . All effective radii are cir- 
cularized, unless noted otherwise. 

2. OBSERVATIONS AND SAMPLE 



1 




Fig. 1. — The WFC3 -ffpi6ow-band image of the galaxy. It is 
well-separated from its nearest neighbors. 



Our study utilises new WFC3/IR i?Fi60W _ band imag- 
ing data taken within the HUDF. These data are part 
of the first of the three ultra-deep pointings which will 
be completed over the next year as part of the HUDF09 
HST Treasury program (G011563). The current WFC3 
imaging consists of 78600 seconds of exposure time in 
the i?Fi60W band, leading to a limiting magnitude of 
28.8. The post-spread function (PSF) FWHM is - 0.16 
arcsec. Details of th e data reduction can be found in 
IBouwens etHI (|2010f ). 

Since the WFC3 data does not cover the complete 
HUDF, m ost of the compac t ma ssive z « 2 galaxies 
from e.g. iDaddi et al.l (|2005h and iCimatti et aTT (|2008f ) 
fall outside of the observed area. From the compact 
z m 2 galaxies inside the WFC3 HUDF image area we 
select the most massive one, located at a = 3 : 32 : 38.12, 
S = —27 : 47 : 49. 63. This galaxy ha s a spectroscopic red- 
shift z = 1.91 (IDaddi et all 120051) , stellar mass M» = 
0.56 x 10 n M Q (jWuvts et all 120081 : Forster Schreiber in 
preparation) , and effective radius r F < 1 kpc in the ob- 
served z band (IDaddi et al.H2005b ICimatti et al.ll2008H . It 
was identified by IDaddi et al.l ((2005) as passively evolv- 
ing based on the BzK criterion. A summary of the 
galaxy's structural parameters is given in Table [TJ An 
image of the galaxy is shown in Figure [TJ It is sufficiently 
separated from its neighbors to prevent contamination of 
its surface brightness profile. 

3. FITTING AND SIZE 

We use the G ALFIT packag e (|Peng et al.ll2002t ) to fit 
two-dimensional Eersic (1963) model profiles convolved 
with the PSF to the observed surface brightness distri- 
bution. This is an essential step in deriving the struc- 
ture of the galaxy, as the FWHM of the PSF of the 
WFC3 images is significant compared to the size of the 



galaxy. We use a PSF extracted from a nearby unsatu- 
rated star and base our masking image on a SExtractor 
(jBertin fc A rnouts 1996) segmentation map. We fit nine 
different models with a fixed Sersic index (n = 1, 2, 9), 
as well as a model where n is a free parameter. 

The effective radii from the Sersic fits range between 
0.42 and 0.48 for Sersic indices varying between n = 1 
and n = 9, with the free-n fit producing a value of 0.43 
kpc (at n = 3.7). The best-fit Sersic profiles are shown 
in Figure [2] Despite the fact that the effective radii are 
rather similar, it is clear that the derived profiles vary 
significantly with n. 

There is no intrinsic reason why galaxies should have 
"perfect" Sersic profiles. Although locally the surface 
brightness profiles of elliptical galaxies are well fitted by 
single Sersic profiles o ver a large range of radii (e.g., 
IKormendv et al.l 120091 ). the situation may be different 
at high redshift. In particular, if elliptical galaxies 
grow by an inside-out process (e.g. . Hopkins et alT 2009a; 
iFeldmann et al.ll2009D . the surface brightness profiles of 
their progenitors may deviate from Sersic profiles. We 
therefore developed a method to derive more robust in- 
tensity profiles, which depend less on the Sersic n param- 
eter used for the fit. 

Our approach is the following: for each Sersic fit, we 
calculate the residual image, which is an image of the 
observed flux minus the PSF-convolved model. We de- 
rive a profile of the residual flux measured along circles 
centered on the galaxy. We add this residual profile to 
the deconvolved model Sersic profile. We note that the 
intrinsic profile is deconvolved for PSF, but the residuals 
are not. This procedure is similar to how the CLEAN 
deconvolution method employ ed in radioastronomy han- 
dles residuals (lHogbomlll974l ). We thus remove or add 
flux at those radii where the model does not adequately 
describe the data, making a first order correction for 
errors caused by the incorrect profile choice. For large 
radii, where (systematic) uncertainties in the sky deter- 
mination become significant, we extrapolate the residual- 
corrected profile by using the uncorrected Sersic profile, 
scaled to the residual-corrected profile at the transition 
radius. These "residual-corrected" profiles are then in- 
tegrated in order to determine the true half-light radius, 
which we refer to as r e ^ econv . The residual-corrected pro- 
files are shown in Figure [2] The structural parameters of 
the best-fitting profile are given in Table [TJ 

It is clear from Figure[2]that the residual-corrected pro- 
files are much less sensitive to the Sersic-n value adopted 
for the initial modeling, especially at radii beyond a few 
kpc. Furthermore, deviations from the Sersic profile are 
taken into account; as we show in Section 4, using the 
residual-corrected profile we can trace the true surface 
brightness profile much more accurately than using sim- 
ple analytical Sersic fits. This is due to the fact that the 
S /N of the faint emission at large radii is so low that the 
fitting procedure ignores it, even though a lot of flux can 
originate there. Thus the stability of the parameters de- 
rived from Sersic fits is no guarantee for correctness. This 
is particularly relevant when the galaxies have complex 
morphologies, such as in the case of a bright, compact 
galaxy surrounded by a faint, extended envelope. 

Uncertainties in r e .d eco „„ and the total 7?-band mag- 
nitude are estimated from the range in values obtained 
from the fixed-n residual-corrected profiles. The errors 
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TABLE 1 

Structural parameters 



Source 


n 


r e (kpc) 


b/a 


M, (lO n M ) 


-^Fieow ( AB ) 


This Letter 
Previous work 


3.7 ±0.38 
4.7 ±0.6* 


0.42 ±0.14 
0.79 ±0.08! 


0.70 
0.74 1 


0.56 2 


22.15 ± 0.067 
22.12 ± 0.03 2 



1 Daddi ct al. (2005), measured in the ZF850LP band 

2 IWuvts et"aLl [I2008T ) 
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Fig. 2. — Our method of correcting the observed surface brightness profile for the effects of the HST WFC3 PSF and incorrect profile 
modeling. In the left panel the best-fitting Sersic models, derived from a two-dimensional fit using a star as the PSF, are shown for different 
values of n. The black curve indicates the free-n fit, with n = 3.7. The profiles show large variations. In the top center panel the observed 
profile is shown. The residual fluxes from the Sersic fits are shown in the bottom panel as a fraction of the observed flux. In the right 
panel the profiles derived using our "residual-correction" method are shown. At large radii, where uncertainties in the sky determination 
become significant, the profile is extrapolated. This is indicated by dashed curves. The residual-corrected profiles are much more robust 
to modeling errors than the uncorrected profiles. The derived effective radius is indicated on the bottom x-axis, the PSF size (HWHM) is 
indicated by the star symbol on the top x-axis. The solid horizontal line in the middle panel indicates the 3<r sky noise level. As can be 
seen, the surface brightness profile can be robustly measured to a surface brightness of 28 mag arcsec -2 . 



given in Table [T] are the rms errors of the best-fit pa- 
rameters from all of the fits, and give an indication of 
the systematic errors due to differences between the ob- 
served surface brightness profile and the Sersic models 
used in the fitting procedure. The uncertainty in n is 
estimated using simulations: we add random sky noise 
to the observed galaxy image. This is repeated several 
times, resulting in a number of images, on each of which 
we perform the fitting procedure described above. The 
uncertainty given in Table [T] is two times the rms error 
of the best-fit parameter from all of the fits. 

Our results are the following: the galaxy is best fit by a 
Sersic profile with n — 3.7. Using the residual-corrected 
profile we find that the effective radius of the galaxy is 
r e ,deconv = 0.050 arcsec, which corresponds to r e ^ econv = 
0.42 kpc. If we fix the Sersic index to a constant value, 
the inferred size does not vary substantially: r e ^ e conv 
varies from 0.31 kpc for n = 9 to 0.51 kpc for n = 1. 
Thus, the deviations from the best-fitting profile are < 
20%. Our size estimate is therefore reasonably robust to 
deviations from the model profile. 

We have investigated the influence of PSF uncertain- 
ties; if we use PSFs extracted from other stars in the field 
we find variations in r e ^econv of < 10%. We have used 
the Tiny Tim software packag^H to investigate the spatial 
dependence of the PSF independently. We find that the 
derived effective radius changes very little with the posi- 
tion of the reference star used, with a maximum of 10% 
in opposite corners of the field. The difference in effec- 

7 http: //www. stsci . edu/sof tware/tlnytim 



tive radius due to the distance between the reference star 
and the galaxy is less than 1%. We therefore conclude 
that PSF errors do not present a significant problem in 
our analysis. 

4. LOW SURFACE BRIGHTNESS SENSITIVITY 

We now determine whether faint extended emission 
would be detected using our data. To this end we con- 
struct several simulated galaxy images which consist of 
two components; a compact component, described by 
a n — 4 Sersic profile with an effective radius roughly 
equal to the observed galaxy (see Table [1]) , and an ex- 
tended component, described by a Sersic profile with ei- 
ther n = 4, r e ?s 3.5 kpc or n = 1, r e ~ 15 kpc. The 
extended component has a flux that is either 10% or 50% 
of the compact component's flux. The compact compo- 
nent's flux is chosen such that the total flux of the two 
components is equal to the observed galaxy's total flux. 
The images are convolved with the PSF, and sky and 
readout noise are added. The images are then fit with 
a single Sersic profile using GALFIT, and a residual- 
corrected profile is constructed. By comparing the half- 
light radii obtained in this way to the intrinsic half-light 
radii we can quantify the sensitivity of our data to low 
surface brightness components. 

The results of our simulated galaxy fits are shown in 
Figure [3] The residual-corrected profiles closely follow 
the intrinsic profiles. The effective radii derived from 
the residual-corrected profiles are very close to the in- 
trinsic effective radii: in three of the cases the difference 
is less than 5%. For the n = 1 extended component with 
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Fig. 3. — Residual-corrected fits to simulated galaxy images (red curves). Top: compact n = 4 profile and extended n . = 4 profile, with 
flux ratios 10 : 1 {left) and 2 : 1 (right). Bottom: compact n = 4 profile and extended n = 1 profile, with the same flux ratios as in the top 
panels. At large radii, where uncertainties in the sky determination become significant, the profile is extrapolated. This is indicated by 
dashed curves. The solid black curves indicate the total intrinsic surface brightness profiles, the dotted black curves indicate the individual 
components that make up these profiles. The best-fit uncorrected Sersic profiles are shown as blue curves. These deviate strongly from 
the true profiles at large radii. The residual-corrected profiles (red curves) follow the intrinsic profiles extremely well, demonstrating that 
our method recovers the intrinsic profiles accurately; the derived effective radii, indicated on the bottom x-axes, are within 10% of the true 
effective radii. The PSF size (HWHM) is indicated by the star symbols on the top x-axes. 



a total flux equal to half of the compact component's 
flux the inferred radius is 10% smaller than the intrinsic 
radius, comparable to the systematic error due to mod- 
eling uncertainties (see Section [3]). We also tested n = 4 
and n = 1 models with effective radii of several kpc for 
the n = 1 extended component: these models are so well 
approximated by Sersic models with higher values of n 
(> 4) that normal Sersic profile fitting immediately re- 
trieves the correct effective radii. 

In conclusion, our method used on these deep data 
is sensitive to a faint extended component down to a 
surface brightness of H as 28 mag arcsec -2 , and using our 
method we retrieve effective radii that are within la of 
the true value. We note that the effective radii obtained 
using the conventional method are, in most cases, very 



close to the intrinsic effective radii. However, the surface 
brightness profiles obtained in this way clearly deviate 
from the intrinsic profiles. 

5. DISCUSSION 

We have found that the galaxy under consideration is 
indeed remarkably small. We have fitted a Sersic model 
to the observed flux distribution, and corrected the pro- 
file for the observed deviations. We have measured the 
galaxy's half light radius: r e ^ eC onv = 0.42 ± 0.14 kpc. 
This result is robust to changes in the imposed Sersic 
profile. As a check of our data's sensitivity to a low 
surface brightness component we have constructed simu- 
lated galaxy images which include a faint extended com- 
ponent. We can reproduce the effective radii to 10% 
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Fig. 4. — Relation s between s i ze and stellar mass (left) and size and rest-frame r-band luminosity (middle) for a sample of low-redshift 
galaxies, taken from fGuo et al.l (2009J). The large symbol with error bars indicates the position of our galaxy. Low-redshift galaxies are 
much larger at similar stellar masses and luminosities. The arrow in the middle plot indicates the change in luminosity due to passive 
evolution to z = 0. The size of the galaxy is smaller than the local equivalents by a factor of 10. Right: compari son of best-fit residual- 
corrected rest- frame V-band surface brightness profile (red curve) to elliptical galaxies in the Virgo cluster, from Kormcndy ct al. (2009) 
(black curves). Virgo galaxies are plotted with M^/Mq > 10 11 . The observed high-z surface brightness profile has been corrected for 
cosmological surface dimming and passive M/L evolution from z = 1.91 to z = (see text). Assuming the galaxy has a mass > 1O 11 M0 
at z = its profile at large radii will evolve very strongly over the next 10 Gyr. The central surface brightness profile, on the other hand, 
shows much less evolution between z 2 and 2 = 0. 



using our technique. 

A possible cause for concern is that the galaxy might 
deviate strongly from a Sersic profile. We have incor- 
porated the residuals in our fit to compensate for such 
errors, and we note that the residuals from our best Sersic 
model fit are quite low (< 10%). This implies that our 
model profile is close to the real profile. This, and the 
fact that varying n has little influence on the derived half- 
light radius, suggests that our results are not strongly 
affected by this source of error. 

Thus, our findings indicate that the small effective ra- 
dius that has been found is not due to oversimplified 
modeling or a lack of S/N, and gives additional evidence 
that a strong evolution in size occurs from z 2 to z = 0. 
It should be noted that our derived effect ive radius is 
1.6 tim es smaller than the radii derived by iDaddi et al.l 
(2005) in the i and z bands. When we repeat our anal- 
ysis on the ACS z-band data we obtain a slightly differ- 
ent value, r ej deconv ~ 0.65 kpc (uncircularized) , closer to 
the deep H-ba,nd i maging, and somewh at smaller than 
the value derived bv IDaddi et al.l (|2005| ) (but consistent 
within the errors). Hence all bands indicate a very small 
size. 

Figure [4] illustrates the difference in size and mass be- 
tween our galaxy and the z = elliptical population; 
plotted in the first two panels are the compact galaxy we 
have studied and a sample of low-redshift central galaxies 
from groups and clusters in the S loan Digital Sky Sur- 
vey, analyzed by IGuo et all (|2009() . The compact z 2 
galaxy lies far off from the z = mass-size relation. The 
middle panel shows the galaxy on the mass-luminosity 
relation. We estimated the luminosity evolution of the 
compact galaxy from z = 1.91 to 2 = in two ways: 
we first used the rest-frame B — I color difference be- 
tween low and high redshift to estimate the difference 
in mass-to-light ratio. Second we used the Fundamental 
Plane to estimate t he ev olution from z — to z — 1 from 
van der Wei et all (|2005[ ) , and used the average evolution 
of the mass-to-light ratios of early-types in the CDFS 
at z = 1 and the z = 1.91 galaxy, both from Forster 
Schreiber et al. (in preparation). The resulting evolu- 
tion is 1.8-2.2 magnitudes. As a result, the galaxy still 



lies off from the size-magnitude relation after correcting 
for evolution. 

In the third panel of Figure |4] we compare the sur- 
face brightness profile of this galaxy to those of ellipti- 
cal galaxies in the Virgo cluster. The profile shown has 
been corrected for cosmological surface brightness dim- 
ming and passive luminosity evolution from z = 1.91 
to z = 0. The total correction is —3.5 + 2 « —1.5 
magnitudes. Even though the galaxy has an average 
density > 100 times larger than the average z = el- 
liptical of the same mass, its surface brightness profile 
in the central kpc is actually rather similar to those of 
the most massive galaxies at z = - the average den- 
sity measured at fixed physical radius is not that differ- 
ent. This is consistent with result s obtained by other au- 
thors (e.g.. iBezanson et al.ll2009t iHopkins et al.l [2009a; 
iFeldmann et al.l 120091: Ivan Dokkum et al.l l2010f T Thus, 
the main difference between z = and this z « 2 galaxy 
is at larger radii where the z ss 2 galaxy has much lower 
surface brightness. Such a result could be explained by 
inside-out growth. 

We also note that there may be significant errors in 
the mass determination of z « 2 compact galaxies, due 
to e.g. incorrect assumptions about the IMF. Changes in 
the low mass end of the IMF affect both the masses of 
the high redshift and low redshift galaxies, and are nearly 
irrelevant. However, changes in the slope of the IMF will 
affect the derived passive evolution between z = 2 and 
z = 0, and will increase or decrease the size evolution. 
Changes in the IMF could thus have important conse- 
quences for evolution. Future deep NIR spectroscopic 
data should provide direct information on the kinemat- 
ics of these ob jec ts and will allow us to con firm their high 
masses (see e.g.. Ivan Dokkum et alll2009l ). 

Finally, it will be interesting to obtain similar deep 
data on other compact massive galaxies, so that their 
profiles can be analyzed to the same surface brightness 
limit. We note that stacking ca n also lead to a grea t 
increase in imaging dept h; e. g.. iCassata et afl (12009). 
Ivan Dokkum et alj pOOcl . and lvan der Wei et all {20081 ) 
stack samples of compact galaxies and obtain very good 
constraints on their average surface brightness profile. 
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ing years many more compact massive galaxies can be 11563 and ERC grant 227749. 
studied on an individual basis. 
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